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Abstract 

In this work, we study some physical aspects of unitary evolution of Bianchi-I model. In particular, we study 
the behavior of the volume and the scale factor as a function of time for the Bianchi-I universe with ultra- 
relativistic fluid (a = 1). The expectation value of volume is shown not to hit any singularity. We elucidate 
on the anisotropic nature of the solution and physically interpret the wavefunction as a superposition of 
collapsing universe and expanding universe mimicking Hartle-Hawking type wavefunction. The same analysis 
has been done for a ^ 1 as well. The work also comments on how to obtain Vilenkin type wavefunction 
representing an expanding universe without that being superposed with a solution representing a collapsing 
universe. We also point at a possible phase transition phenomenon for a ^ 1 case with existence of a critical 
parameter. 
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1. Introduction 

Quantum Cosmology, an amalgamation of Quantum theory with Cosmology, deals with quantizing a 
cosmological model as a whole. The quantization scheme builds upon the famous Wheeler-DeWitt Equation 
[1, 2] paving the way to have a hold on the physics of the beginning of the universe. However, we have 
several conceptual and technical problems associated with the formulation. One of the important problems 
involves suitably choosing a time like parameter since in a relativistic theory, time itself is a coordinate 
[3, 4, 5]. In another study [6], a scheme has been proposed where time plays no role whatsoever. Apart from 
this, there are conceptual problems regarding the interpretation of the wave function, implementation of 
proper boundary conditions, and several others. The various aspect of these subtleties have been reviewed 
along with the development of the subject in [7, 8, 9]. 

The issue of conjuring up a time parameter can be resolved if matter sector comprises a perfect fluid, a 
scenario, which enables us to exploit a formalism, relating the four velocity of the fluid with thermodynamic 
parameters, developed by Schutz [10, 11]. It so happens that after some clever canonical transformations one 
can show that the canonical momentum associated with fluid variable appears as a linear term in Hamilto¬ 
nian and thereby mimics Schrodinger equation. Hence the variable canonically conjugate to this momentum 
can be treated as the time parameter. Lapchinskii and Rubakov[12] utilized this strategy to identify a well 
behaved time parameter. This method has found a renewed and very successful application by Alvarenga 
and Lemos[13] and has subsequently been used, for the quantization of an isotropic cosmological model, in 
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the works of Batista et al[14], Alvarenga et al[15] and Vakili[16, 17]. The method also finds an application 
in quantizing anisotropic models in the work done by Alvarenga et al[18], Majumder and Banerjee[19], Pal 
and Banerjee[20]. 

One of the pathology that came up while employing the method to anisotropic universe is the non uni- 
tarity associated with evolution operator, which is a result of having a non self-adjoint Hamiltonian. This 
pathology was present in the Bianchi I model investigated by Alvarenga et. al[18]. It should be mentioned 
that without matter sector, we do not have any physical identification of time, hence the non unitarity 
associated with the model might escape notice [21, 22] while the use of fluid evolution enables us to study 
the issue critically. 

Majumder and Banerjee [19] showed that a suitable operator ordering can make the norm of the wave 
function, corresponding to an anisotropic model, time independent in the asymptotic limit. Although as a 
whole the model still suffered from non unitarity, it was a good enough indication that the problem of non 
unitarity has to do with operator ordering rather than being a generic feature. The confirmation comes in 
a later investigation[20], where it has been shown, contrary to the folklore, that the non unitarity is indeed 
not a generic feature of anisotropic cosmologies. With a proper choice of operator ordering, unitarity can 
be restored to the model. Initially it was shown for Bianchi-I model, later the work has been extended to 
encompass the Bianchi-V and Bianchi-IX universe as well[23]. 

This work sheds light on some physical aspects of “Unitarity Restored Bianchi-I Model” building upon 
the work in [20, 23]. We elucidate on the nature of wavefunction along with a critical study on evolution of 
scale factors, volume, manifest anisotropy of the solutions to the model. The physical insight strengthens the 
idea presented earlier in [20, 23] and leads to a better understanding of the physics of a quantized anisotropic 
model. 

The paper is organized as follows. In section two, Schutz formalism is briefly described along with 
quantization of Bianchi-I model and a discussion on the evolution of expectation value of volume, the scale 
factors associated with the model. In Section three, we take up the a = 1 case in detail and show that it is 
possible to have a quantized unitary Bianchi-I model with an explicit anisotropy. In section four, we discuss 
the a ^ 1 case, elucidate on the nature of the wavefunction and the eigenspectra of energy. We also point 
at a possible phase transition with existence of a critical parameter. The fifth and final section discusses 
the results obtained in the work and its implications. 


2. Quantization of Bianchi-I Cosmology 


The general action for gravity is given by 

A = [ d*x^R + 2 [ VhhabK‘^’‘+ f d^x^P, (1) 

JM JdM Jm 

where K°‘^ is the extrinsic curvature, and /i®*” is the induced metric over the boundary dM of the 4 dimen¬ 
sional space-time manifold M. The units are so chosen that IGttG = 1. 

The metric for Bianchi-I is given by 

ds^ = n^df - [a^{t)dx^ + h^{t)dy'^ + c^{t)dz^] , ( 2 ) 


where n{t) is the lapse function and a, b, c are functions of the cosmic time t. With this metric the gravity 
sector of (1) takes the form 

2 


Ag = I dt 


I “‘Hi 


-( abc + bca + cab 

n 


( 3 ) 
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With a set of new variables 


a{t) = e^o+(i++V3fi-^ 

b{t) = , (5) 

c(t) = (6) 

the Lagrangian density of the gravity sector becomes 

C, = -6^[p^o-P--Pl)- (7) 

The corresponding Hamiltonian density can be written as 

=-nexp(-3^o) (Po-P+-P-)| • (8) 


Now the Schutz’s formalism [10, 11] is employed to identify a time parameter out of the matter sector, which 
comprises an ideal fluid given by an equation of state P = ap where a is a constant (with a < 1), P and p 
is the pressure and density of the fluid respectively. Standard thermodynamical consideration leads to the 
following action for the matter sector, 



dtCf 


= V 


dt 




(1 + a) 


1+^ 



(9) 


where 0, e and S is defined as in the reference [20, 23]. We define canonical momentum to be Pe = and 
PS = and the corresponding Hamiltonian density comes out to be 


Hf = 

With the following canonical transformation, 

T = -psexp{-S)p~°‘~^, 

PT = p“+^exp(S'), 

' _ , f , PS 

6 — 6 + (q + 1) -, 

Pe 

P'e = 


( 10 ) 


( 11 ) 

( 12 ) 

(13) 

(14) 


which maintains the canonical structure by preserving the Poisson commutation relation[20], the Hamilto¬ 
nian density for the fluid sector can be cast into following form 

Hf = (15) 

Combining gravity and fluid sector, the super Hamiltonian (to be precise, super Hamiltonian density) 
becomes 

H = (pI -pI- pi - (16) 

and variation with respect to n yields Hamiltonian constraint 

n=- = -{Hg + Hf) = {). (17) 

nn 
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We now promote the variation of super Hamiltonian with respect to n i.e % to an operator and postulate 
commutator relation amongst the quantum operators as usual i.e we have: 


for i = 0, +, and 

This mapping is equivalent to postulating 


Pj -ihdj3., 

Pt — iHOt- 

the fundamental commutation relations: 


(18) 

(19) 


[l3j,Pk] = (20) 

In all subsequent discussion, we employ the choice of natural unit i.e H = c = 1 and rewrite the equation 
(17) in its operator form as T-Lip = 0, which can be recast, using the operator ordering as prescribed in [20], 
in following form: 




a/3o 


(“-l)do_ 




3(a-l)/3o 


52 

dpl 


„3(a-l)/3o 


92 

dpi 




( 21 ) 


Stiff Fluid P = p 

Now for ultra relativistic fluid we have a = 1 and the equation (21) reduces to 

q2 q2 q2 

% 


dpl dPl dpi 


d 

ih = 24:1-—lb. 


( 22 ) 


Naively, (22) seems to be non self-adjoint. But we can calculate deficiency index (n±) of the operator, 
which turns out to be 0, hence the operator is self-adjoint and the corresponding evolution is guaranteed 
to be unitary. The general solution to eq. (22) has already been described in [20]. In order to extract the 
physical informations, here we consider the solution corresponding to an eigenstate of both p_|_ and p- with 
eigenvalues = k- = 0, (we will see, what it does mean physically to have = 0 as we go along) which 
is given by 

if = ^ 23 ) 

2 

where energy, E is given hy E = — Now using linearity of the governing differential equation we obtain 
the following wave packet 

r {u: + 12u:„p 

* ipau), 




which upon integration and normalization, can be recast as 




-36<.= / 

= exp -- 


l-ii 


l-*6 


(24) 


(25) 


with a time dependent probability density function, given by 


= 


V2 


^(1 + ie) 


: exp 


-2 iPo - iOoTf 

, 1 + il 


which, on integration over the conhguration space, yields a time independent value: 

J '4/*4/dPodP+dP- “ 4 y dP-dP+ = l.d (0) S (0). 


(26) 


(27) 


It deserves mention that presence of <5(0) indicates that the wavefunction belongs to the rigged Hilbert space 
just like the eigenfunction of momentum in elementary quantum mechanics. 
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The definition of norm used here (27) is consistent with the definition of norm used for a ^ 1 case, in 
the earlier work [ 20 ], in the sense, by properly taking limit a —>■ 1 the norm used there reduce to the one, 
used above. For a 7 ^ 1, we have, for /3o sector, [20] 

IIV’IU/i = J (28) 

whereas for a = 1 , the norm is defined as 

llV-lla^i = / = lim 3 (29) 

We note, the limit is well-defined since, (28) has a factor of (1 — a). This also suggests that for a 7 ^ 1, 
we could have defined norm without factor, because physics is independent of such redefinition of a 

norm by a numerical factor. In that scenario, we don’t even need to divide by the factor in (2) to 

get the definition of norm for a = 1 case from a 7 ^ 1 . 

Having defined norm in proper way, now, for the wavefunction (25), we can evaluate the expectation 
value of /3o, which behaves as follows: 

(^ 0 ) = tooT. (30) 

while the expectation value of volume is given by: 

(F) = (31) 

If we set Wo = 0, the volume expansion can be depicted via figure 1. Although (/3o) is 0 in this case (i.e 
with Wo = 0), still we have expansion due to quantum fluctuation in the variable Pq. Another important 
aspect is that the volume assumes a finite value at T=0 and thereby avoids hitting the singularity as evident 
from figure 1. It deserves mention that by setting wo = 0 we do not loose any qualitative feature. As long 
as Wo > 0, we do get an expansion. Only if wo becomes negative, we would get a contracting phase until 
T = 48|wo|, after which, the universe will start expanding. 


Spatial Volume Expansion for <^>=0 
V 



Figure 1; Behaviour of Volume with time for (do) = 0 


Looking at the behavior of the wavefunction with time, as shown in figure 2, we observe that the peak of 
the Gaussian shifts toward right and spreads out for T > 0. The spreading out happens because within the 
superposed solution there are different w modes having different “velocity”, which ultimately leads to dis¬ 
persion. The term “velocity” corresponds to the rate of change of volume since the conjugate momentum to 
/3q is actually proportional to the time derivative of volume factor. Thus, different w modes may be assigned 
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Gaussian WavePacketiProbability Density for </3o>=T 



Figure 2: The Probability Distribution as a function of 00 at different time T 


the physical interpretation of having different values for the rate of change of volume. This scenario can be 
compared to a Gaussian wavepacket of a free particle dispersing, in quantum mechanics. We note that the 
choice of Gaussian prefactor is based on the fact that it leads to a behaviour where the universe does not 
hit a singularity at T=0, we get expansion in the volume of the universe and we can nicely interpret the 
behaviour of the wavefunction. 


The wavefunction, given by (25), is basically Hartle-Hawking type, since we have superposition of both 
positive and negative ui modes. We can think of positive uj modes to describe an expanding universe whereas 
the negative modes can be thought of depicting a collapsing universe. Since these lu modes are weighted by 
Gaussian pre-factors, most of them gets suppressed and effectively we have an expanding universe. 

As a side note, we observe dt is also a solution to isotropic quantum FRW model with a spatially flat 
metric since is independent of /3+, /3_. Moreover, if we restrict to isotropic FRW model, we do not need 
to integrate over /3+, /?_, hence the delta functions would not appear in the norm (27) and the wavefunction 
(25) perfectly belongs to the Hilbert space. 


The fact that does satisfy Hamiltonian constraint for quantum FRW universe, can be understood from 
inspecting /3o,± in terms of scale factor. We note that 


: -ln{abc) = -InV, 

(32) 

1 f o,b\ 

(33) 

, 

(») ^ 

(34) 


Hence, /3o actually dictates the volume of the universe while the anisotropy is dictated by the variables j3± 
and variables conjugate to them. If we have isotropic solution, the scale factors will have similar behaviour 
and we expect k± = 0 and (/3±) = 0. Here we note that the p± sector admits plane wave like solution i.e 
gi(fc+/3++fc_/3_) _ Q have (/3±) = 0 where the following limiting procedure is implicitly 

assumed 




( 35 ) 
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. Hence, the solution is indeed isotropic in nature. The fact that the anisotropy is dictated by j3± and p± 
is evident if we look at the scalar found out of shear tensor 


1 

12 


( ^11 

V5ii 



1 (P\ 

12e6/5o 1^18 



( ^ 
V522 



/to _ gii V 

Vs'33 511/ 


(36) 


Now one way to make the model manifestly anisotropic is to take (/3±) 7 ^ 0 or to make the momentum 
conjugate to j3± nonzero. It deserves mention in this connection that the choice of /c± = 0 may seem to 
be restricted, but in the next section, we will relax the condition and make k+ 7 ^ 0. To be precise, we will 
construct solution with nonzero (/3+) and In fact, it is quite trivial to extend the analysis to k± 7 ^ 0, 
we will get similar Gaussian solutions for /3± sector as well. We have not done it since the essential physical 
features are conveyed with the restricted condition on k ±. This kind of choice has also been made earlier in 
the work of Alvarenga et. al. [18]. 


3. Explicit Anisotropy 

In the last section we constructed an isotropic solution. Hence, one can argue that the unitarity is 
obtained at the expense of anisotropy. In this section we present a solution, which is manifestly anisotropic. 
We consider eigenstate of p- with eigenvalue fc_ = 0 but let lo and to vary along with energy eigenvalue 
given by E = — ^ {uP' — fc^). From the general solution ip = e*G/3o+fc+/3+-BT)^ construct the following 
superposed solution : 

T f f ('°+-12'=+o)^ /o^\ 

vl/=//e ^ e ^ 'ipdujdk+, (37) 

which, upon integration and normalization, yields 


/ 2 p-36(i^^+k^c) (/3+-.6fc+o)2 

= W- c i-'T e 


1 + 


36 


with a time dependent probability density function, given by 


vE'*^' = 


-(1 + S) 


exp 


-2 [(do - + (/3+ - k+oT)^] 


1 + ^ 
^ 36 


(38) 


(39) 


which leads to a time independent norm, as evident from 


d/3od/3+d/3-^*^ = 1(5(0). 


(40) 


The (5(0) is there as we have a plane wave like solution for k_ sector. Now, (38) has the property 

(do) = uJoT, (/?+) = k+oT (41) 

(d-) = 0 (42) 

We have thus constructed an anisotropic solution explicitly since we have nonzero value for (d+) and 
superposition of nonzero fc_|_ modes. Having (d-) = 0 means that the scale factor a and b behaves similarly 
while non-zero expectation of d-i- means scale factor c(t) behaves differently, as evident from equation (33). 
The signature of anisotropy is also carried by expectation value of which satisfies. 


(d^) = 


-6/3o^ 


12 


18 8 


(e-6/3o) 

216 


PO 


(43) 



















While constructing the wavepacket (38) , we have not put the constraint on the energy spectrum that 
it has to be negative on physical ground because we have Hamiltonian constraint (17) and we expect the 
energy of fluid to be positive definite[20]. Now negative energy yields the constraint This, in turn, 

implies integration over a restricted interval in (A:+, w) space while obtaining the superposed solution, which 
would be difficult to inspect analytically. Nonetheless, the norm for this case remains finite. 


and is given by: 


4. A more general fluid with a ^ 1 

The solution to the (21) for a 1 has already been described in the reference 

3{1-q:)/3o 

where (pa^b are given in variable x = e 5 by 

Ux) = VX (Ax) + (Ax)] , (45) 

Mx) = Vx (Xx) + (Ax)] , (46) 

respectively for cr > |; and cr < |, where g = —la = ^a — \ with a = ^ so that 4>a,b satisfy. 


(P(f) 

dx^ 



-E'p. 


(47) 


In both cases, the spectrum is given by E' = E = — A^. The self-adjoint extension guarantees that 

^ takes a value so as to conserve probability and makes the model unitary [20]. For a 1, it is technically 
difficult to find the expectation value of Pq due to mathematical complexity of the integral involved 


iPo) = 


J dxln (x) |(/>|^ 


3(1-0) J dxlpr^ 


(48) 


where (p is given by (pa or (pb for cr > ^ or u < ^ respectively. Nonetheless, as we are dealing with energy 
eigenstate, this expectation value ought to be time independent trivially. 


We now try to extract some physical information even without the integrated version. Asymptotically, 
the solutions (45), (46) become a linear combination of incoming and outgoing waves. Since, the eigenvalue 
of (momentum conjugate to the newly defined variable x) is proportional to rate of change of volume 
factor, , incoming wave can be thought of describing a collapsing universe with negative rate of change 
of volume whereas the outgoing wave describes an expanding universe. Hence, the solution (44), being a 
superposition of expanding and collapsing universe, mimics Hartle-Hawking wavefunction [24]. The prob¬ 
ability conservation can now be physically thought of as the probability of collapsing being same as the 
probability of birth of a universe. 


We can have Vilenkin type solution as well, all we need to do is to keep the outgoing solution while 
discarding the incoming wave solution so that the wavefunction becomes 


(t>a{x) = 'Jx 


BII[^^{\X) 


Mx) = M 


bbA)(\x) 


(49) 

(50) 















Previously, the boundary condition has imposed restriction on the ratio of A and B. Now since, we have 
A = 0, the boundary condition enforces constraint on the energy spectra. It can be shown for a < 1/4, the 
condition coming from self-adjoint extension [25, 20], given by 


TT 



-I- lain 



-I- arg[T (1 -|- a)] = nir, 


(51) 


(where xq characterizes the self-adjoint extension, elucidated in [25]) is satisfied if 


a = 


Xo 

II 

IT 

O 

(52) 

(^) 

= 1, 

(53) 

a (1-0) 

= 2n, 

(54) 

4 (fc^ -1- 

1 4n^ 

(55) 

9(l-a)^ 

4 (1-0)^' 


The condition (53) implies existence of unique energy eigenstate with energy eigenvalue, given by 


E' = 


4 


(56) 


Similar analysis reveals that for a > ^, 6 must be equal to 1 i.e Xq = which yields a unique energy 
eigenstate with E' = — for cr = ^ (where ln{"f) is Euler’s constant) while for cr > j, we have a bunch 

of energy eigenvalues E’ = satisfying 


gin 



- arg[r (1 -F ig)] 


TITT. 


(57) 


Unphysical positive energy states 

In [20], it has been argued that observed universe with positive energy implies that the fluid sector has 
negative energy (because of the Hamiltonian constraint (17)), which is not physical. Here we will argue from 
the point of view of gravity why the positive energy eigenvalue does not make sense physically. 

We have the following solution with positive eigenvalue for cr > j : 


, /2sinh(7ro) , 

W— ^^VxK.9 (Ax). 


(58) 


where g = | (the wavefunction is shown in figure 3 below) with energy eigenvalue given by E' = 

The solution (58) gets exponentially damped for high y (i.e high /?o) value. It implies that we have 
exponentially damped probability for having larger /3q i.e larger volume which is not physical for an observed 
expanding universe. The same conclusion can be drawn even if we do not have the explicit form of the 
wavefunction but we look at governing differential equation (47) in asymptotic limit and realise that for 
arbitrary large y , it becomes 

- ■ (59) 




which does imply exponential decay in y. 


This comes as a sanity check for our theory. Positive eigenvalue for the gravity sector implies negative 
fluid energy which indicates a ‘trapped’ universe, where the volume has a very low probability to get larger 
than some specihed value which depends on the energy. It deserves mention that cr < j does not admit any 
positive energy eigenvalue [25] and this fact hints at we have two different regime of solution around the 
critical point cr = f [26, 27]. 
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Positive Energy Eigenvalue 


<Pix) 



Figure 3: Positive Energy Eigenstate with E' = 1, g = G 


5. Conclusion 

The work reinforces the idea presented in [20, 23] that the alleged non unitarity in a quantum anisotropic 
cosmological model has nothing to do with hyperbolic nature of kinetic term, rather it is intimately linked 
to co-ordinate choice, which in turn affects the operator ordering. Although we do not know the physical 
reason behind why we do have to choose this ordering, yet, with proper ordering, we restore the unitarity and 
thereby make the model worthy of being investigated critically. This, in turn, enables us to have physical 
insight and predictions from this model, which is cured from alleged non unitarity. 

It has been shown in this work that the universe does not hit any singularity at T = 0. The work sheds 
light on the nature of variables j3± and p ±, which encode the anisotropic nature of the solution and thereby 
exhibit that the solution obtained is indeed anisotropic in nature. Furthermore, we physically interpret the 
wavefunction in terms of collapsing and expanding universe and thereby draw a possible parallelism between 
our wavefunction with Hartle-Hawking type & Vilenkin type wavefunctions. 

The solution for a 7 ^ 1, given in reference [20], is shown to be a Hartle-Hawking type wavefunction. 
In this work, we extend our analysis and construct a Vilenkin type wavefunction for which the self-adjoint 
extension puts constraint on energy eigenspectra. For a > j, we have a bunch of eigenvalues while for 
cr < j case, we have a unique solution. The study of this transition with possible breaking of conformal 
symmetry when the system passes through the critical point cr = | is beyond the scope of this work yet it 
deserves mention that similar kind of phase transition associated with inverse square potential has already 
been known[26] and mimics the BKT phase transition [27]. Moreover, the work gives another perspective 
to the physical impossibility of having positive energy eigenstate of Hamiltonian of gravity sector in terms 
of cosmic evolution. 

We also conclusively show that the anisotropy is dictated by the variables /3± whereas the ordering am¬ 
biguity is in the variable /3q, which dictates the evolution of volume. Hence, it lends support to the idea 
presented in [20, 23] that the pathology of non unitarity is not a generic feature of anisotropic models, it 
may crop up even in isotropic case [23], the problem rather lies in choosing the appropriate co-ordinate and 
operator ordering for the variable which describes the evolution of volume. 
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